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A pseudo-Jahn-Teller model describing central atom distortions is proposed for endohedral 
fullerenes of the form AQCeo where A is either a rare gas or a metal atom. A critical (dimen- 
sionless) coupling gc is found, below which the symmetric configuration is stable and above which 
inversion symmetry is broken. Vibronic parameters are given for selected endohedral fullerenes. 
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Shortly after it was established that fullerenes are 
closed-cage structures of carbon, investigators speculated 
that atoms and clusters might be coaxed to occupy the 
space inside the cage. The existence of such endohedral 
fullerenes is now well-established by experiment [Q . Such 
endohedral doping of Cgo shows some promise as a tech- 
nique to modify fuUerene properties such as electric po- 
larizibility, optical and infrared absorption, and magnetic 
susceptibility. 

Several numerical studies have concluded that for 
the case of metal endohedral dopants, such as NaQCgo, 
LiQCgo, and CaQCgo, the metal atom transfers its va- 
lence electrons to the Ceo cage, and establishes an off- 
center equilibrium position, resulting in a net electric 
dipole moment for the molecule. For the case of rare 
gas dopants, however, the central atom sits at the cage 
center. 

A classical model was used by Erwin 0| to explain the 
two disparate cases; for metal dopants, the charge trans- 
ferred to the cage occupies a highly delocalized tin or- 
bital. Erwin treats the resulting negatively charged cage 
as a perfectly conducting sphere. The total energy of the 
metal cation and the conducting sphere is minimized by 
placing the cation off-center where electrostatic interac- 
tion between the cation and the cage can be made more 
attractive as the charge on the cage is polarized. The 
distortion is stabilized in Erwin's model by including a 
short-ranged repulsive interaction of Lennard- Jones form 
for the electron overlap between the cation and cage elec- 
trons. 

Such a classical model would predict that for the case 
of metal dopants, inversion symmetry is broken as the 
central atom is stabilized in an off-center equilibrium po- 
sition; while for closed-shell rare gas atoms, there will 
be no charge transfer and, consequently, no central atom 
distortion. However, this classical model suffers from two 
serious problems. First, for the case of singly-ionized al- 
kali atoms, such a classical model would again predict 
no charge transfer, as such an atom is isoelectronic to a 
rare gas atom. Thus, charged metallofuUerenes such as 
Li"'"@C6o would have the cation symmetrically positioned 
according to the classical model. This, however, contra- 



dicts electronic structure calculations iQi^- Second, for 
the distortive case of metallofuUerenes, the net dipole 
moment calculated from the classical model is zero, as 
the contributions from the cation and the polarized cage 
exactly cancel. 

The purpose of this work is two-fold: (1) to remedy 
the shortcomings of the classical model, and (2) to give 
insight into the physical mechanisms underlying previous 
numerical studies. We present here a model for the dis- 
tortion mechanism where the electronic states are treated 
quantum mechanically. The propensity of the central 
atom for distortion can be understood as arising from 
the pseudo-Jahn-Teller (PJT) effect, where excited elec- 
tronic states are mixed with highest occupied molecular 
orbital (HOMO) through a symmetry-breaking distor- 
tion. We note that PJT models have been successfully 
invoked in other contexts to describe, for example, dis- 
tortions in the case of the F-center |^ and in oxygen and 
nitrogen impurities in silicon ||^. 

We give analytical expressions for the magnitude of the 
distortion in this PJT model, and we find that there is 
a critical vibronic coupling gc] for values oi g > gc, the 
dopant moves off-center, and for g < g^, the dopant sits 
at the cage center. We also find that for linear coupling 
the adiabatic potential surface is spherically symmetric, 
in good agreement with the numerical calculations , 
while the symmetry of the perturbative warping of the 
adiabatic potential can be explained with the addition of 
higher order terms. 

We detail a model for the PJT coupling that involves 
a displacement of the central ion in endohedrally doped 
Geo- We start by describing the electronic structure of 
the Geo molecule. The valence electrons that are of in- 
terest are the radially-directed p electrons on the car- 
bon atoms that take part in tt bonding. There is one 
such electron on each carbon atom, the remaining va- 
lence electrons taking part in a bonding and having a 
much lower energy. From the atomic p orbitals on each 
carbon, we construct molecular orbitals consistent with 
the Ih symmetry of the Ceo molecule. However, since 
Geo is nearly spherical , it is possible to associate the 
different manifolds of molecular orbitals with angular mo- 
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mentum states L. The Ceo cage is electrically neutral; 
hence, the HOMO is the filled hu level coming from the 
L = 5 molecular orbitals The lowest unoccupied 

molecular orbital (LUMO) is a tiu manifold, also coming 
from L — b. Slightly higher in energy however is a tig 
manifold, from L = 6. It is this excited state that plays 
a crucial role in our PJT coupling model. 

We now consider the Jahn- Teller (JT) interaction in- 
volving the displacement of the central ion from its on- 
center position. For the case of dopants isoelectronic to 
the rare gases, the HOMO remains the filled hu manifold. 
Such a state is non-degenerate, and consequently there is 
no pure JT interaction. There is, however, a PJT interac- 
tion involving the excited tig manifold. The tig manifold 
is the one of interest since the displacement modes of the 
central ion Q span the irreducible representation ti„, and 
tig is contained in the direct product /i„ ® 

The Hamiltonian is taken to be 



H ~Ho+ 'Helas + HjT 



(1) 



TLo is the electronic Hamiltonian and gives the electronic 
energy of the h^. and tig manifolds for the fully symmetric 
(1,1 ) system 



Wo=0 \K){K\+^ \tlg){tlg\, 



(2) 



where A is the energy difference between the tig and 
hu manifolds, A = ^{tig) — E(/iu). Tieias is the elastic 
energy Hamiltonian for a classical displacement of the 
central ion. Assuming the harmonic approximation. 



"Helas = 2 {Qx +Qy + Q z) I: 
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and V gives the strength of the PJT coupling. We note 
that h^^ = a/3/8 he — \/5/8 /le and hx2_y2 = ^5/8 hg + 
■\/3/8 hf^ , in which hg and /le are the states employed in 
Ref. p. 

We proceed by diagonalizing the Hamiltonian in Eq. 
By way of example, we consider the simple case in which 
Qz = Q and Qx = Qy = 0. That is, we consider the 
particular case in which the central ion displaces along 
one of the C2 symmetry axes. In this instance, the elec- 
tronic Hamiltonian matrix reduces to three 2x2 blocks 
and two 1x1 blocks. The h^2_y2 and the /i^ states are 
unaffected by this particular distortion; and so, their en- 
ergies remain unchanged. One 2x2 block, which mixes 
the states h^^ and tiz, has the form 



cVQ 
cVQ A 



(6) 



where c = 2. 

The two remaining 2x2 blocks, mixing the pairs of 
states {h^,tix} or {h^,tiy}, also have the form given in 
Eq. I with c = V3. 

The energy shift for a distortion Q is found by sum- 
ming over the ten lowest energy states, including spin 
degeneracy. Hence, we find 



where I is the electronic identity operator. Actually, Q 
includes ti„-symmetric distortions of the Cgo cage as well 
as a displacement of the central ion. But since the cage 
is much more massive than the dopant atoms considered 
(large cage mass approximation) , the collective displace- 
ment Q is well approximated by dopant displacement 
alone. Lastly, we have the H.jt that describes the PJT 
coupling between the hu and tig manifolds resulting from 
a displacement of the central ion from its on-center po- 
sition, {hu ® tig) (g) tiu coupling. These matrix elements 
are determined from the coupling coefficients given by 
Fowler and Ceulemans (Table 5 of Ref. |12|) where the co- 
ordinate axes are chosen to be coincident with the three 
orthogonal C2 rotation axes. In our eight dimensional 
electronic space, we have 
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where 
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[(t)^ + 3 Q^] Wc then minimize E with respect 
to Q. We find that Qo = 0, corresponding to no dis- 
placement of the central ion, gives the minimum energy 
configuration when JT coupling is weak. More precisely, 
if 



fcA - 40' 



(8) 



there is no displacement of the central ion. On the 
other hand, if the JT coupling is strong enough, g > j^, 
the minimum energy configuration is found for Qo 7^ 0. 
That is, the total energy of the endohedrally doped Cgo 
molecule can be reduced if the central ion displaces from 
its on-center position. The minimum energy configura- 
tion is found when Qo = V Qo / ^ satisfies the equation 



(9) 
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We now consider the most general displacement of the 
central ion. Here, we parametrize Qx, Qy, and Qz by 
Qx — QsmO cos(/), Qy = QsmO sin0, and — QcosO. 
We substitute this into the matrix for the full Hamilto- 
nian. The corresponding characteristic equation, whose 
roots yield the energy eigenvalues, is found to be indepen- 
dent of 9 and </>. Thus, the one-electron energy eigenval- 
ues are independent of the direction of the displacement 
of the central ion, and the adiabatic energy surface is 
spherical. 

Such an accidental SO (3) symmetry has been known 
to occur in other Jahn- Teller systems with icosahedral 
symmetry [|l^. The adiabatic potential E as & func- 
tion of central atom displacement is plotted in Fig. |l| 
for g > gc- When the kinetic energy associated with 
the displacement Q is included in Eq. ^ low frequency 
excitations corresponding to a rotation of the distortion 
about the symmetric configuration are found, together 
with higher frequency radial excitations. Such excita- 
tions have been previously discussed elsewhere It 
is important to note that there is no geometric phase in 
this model. Thus, this model establishes the quantization 
condition for the rotational wavefunctions previously as- 
sumed 1,1. 




FIG. 1. Adiabatic potential iJ vs. magnitude of tiu distor- 
tion Q for Na+QCeo- A second curve (dashed) with reduced 
vibronic coupling V [g = 0.02) is shown where the symmetric 
configuration is stable. 

For the case of metallic dopants, the HOMO becomes 
the tiu manifold [Q. Here, such a partially filled state is 
subject to JT distortions. In fact, it has been proposed 
that such coupling is responsible for superconductivity in 
the alkali-doped fuUeridcs However, the JT distor- 
tion is of hg symmetry, and thus does not contribute to 
a dipole moment. 

In addition to the ii„ ® hg JT coupling, there is PJT 
coupling between the ti„ and tig electronic states by a 



tiu distortion, {tiu © tig) ® ti„ coupling. We take the 
PJT coupling to be 
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Diagonalizing the electronic Hamiltonian yields the 
following eigenvalues: 0, A', and E± = ^ 



1/2 



(^)^ -t- V'^Q^ where A' is the ii„-tig energy gap. 

Again, the eigenvalues are invariant under rotations 
For alkali metal endohedral fuUerenes with a 
ground state, we find that the energy shift for a distortion 
Q is 
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By minimizing the total energy, we conclude that there 
is no distortion unless g > 

A non-vanishing electric dipole moment results from 
the PJT distortion. Because of SO (3) invariance, we con- 
sider without loss of generality the case of a simple axial 
distortion {Qz ^ 0) of an endohedral dopant isoelectronic 
to a rare gas atom. We find that a spontaneous dipole 
moment is formed for a non- vanishing distortion Qz-, in 
contrast to the classical model of Erwin which assumes 
perfectly metallic screening. Using the vibronic eigen- 
states of Eq. |^ together with the Clebsch-Gordan coeffi- 
cients obtained from Ref. |]l2|, we find an electric dipole 
moment in the direction of the central atom distortion of 
the form 



V = 8po 
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where po is a reduced matrix element, a 

Vl + 16Qg-l 
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That numerical calcu- 



lations on light rare gas endohedrals find the equilibrium 
of the rare gas dopant at the cage center indicates that 
g for these systems is less than In the case of posi- 
tively charged alkali endohedrals such as Na+QCgo and 
Li+@C6o, previous calculations ||^,^ have found an off- 
center equilibrium position for the alkali ion, indicating 
that g for these endohedrals exceeds 

Parameters for selected endohedrals are given in Ta- 
ble I For Na+@C6o, we take A = 2.8 eV, Qq = O.TA, 
and E = -0.12 cV |]. From Eqs. and |, we fit the re- 
sults of Ref. |] with V = 0.69 eV/A and k = 5.8 eV/A^, 
giving g — 0.03. 

Assuming the primary coupling is between tiu and tig, 
we estimate the relevant parameters for NaQCgo on the 
basis of previous calculations. We take A' — 1.0 eV, 
Qo = 1.5A, and E = -1.4 eV [|. We fit the results with 
V = 2.47 eV/A and k = 1.63 eV/A^, giving g = 3.73. 
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Ca@C6o can be treated within the same framework. 
For CaQCeo, the electronic ground state for Ih symme- 
try is ^Tig. PJT coupHng with the ^Ti„ state gives rise 
to a distortion that breaks inversion symmetry. The cal- 
culations of Ref. Q can be fit with y = 1.65 eV/A and 
k ^ 3.62 eV/A^, giving g = 0.58. 

The SO (3) invariance of the adiabatic potential within 
the linear PJT coupling will be broken when higher order 
couplings are taken into account. The amount of the 
warping of the adiabatic potential has been found Q to 
be negligible when the central atom dopant is either a 
rare gas or a positively charged alkali atom. Thus, the 
linear PJT model well describes the previous numerical 
calculations for these cases. 

For the case of CaQCeo, it has been found Q that the 
molecular symmetry is broken down to C^y. While the 
addition oi^Tig(^hg linear JT coupling for CaQCgo does 
not break S0(3) invariance, Ds^ warping ||l^ of the adi- 
abatic potential surface does result from quadratic JT 
coupling involving distortions of hg symmetry |l^ , lead- 
ing to a molecular symmetry of Cs^ . The softest modes 
of Ceo are of hg symmetry; consequently, there is a small 
energy cost for such distortions. For such a quadratic 
coupling to the low frequency hg vibrations, the direc- 
tions of three-fold symmetry can contain relative max- 
ima, while the directions of two-fold symmetry contain 
saddle points. For weak coupling, such a warping will be 
small compared to the PJT stabilization energy. 

We conclude that PJT coupling provides a unifying 
framework for understanding the spontaneous polariza- 
tion of cndohedral fuUerenes. The parameters of our PJT 
model were obtained by fitting to previous numerical cal- 
culations. We find that the structural distortion result- 
ing from endohedral doping is yet another example illus- 
trating the importance of vibronic effects in shaping the 
properties of fullerene systems. 
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TABLE L Vibronic parameters for selected endohedrals. 



Dopant 


A or A' (eV) 


Qo (A) 


V (cV/A) 


k (eV/A^) 
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5.78 


0.03 


Na 


l.O'' 


1.5" 


2.47 


1.63 


3.73 


Ca 


1.31*^ 


0.83*= 


1.65 


3.62 


0.58 
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